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We find that the Debye-Huckel nonideality of dilute aqueous electrolytes is sufficient to drive volume phase 
transitions and criticality, even in the absence of a self-attracting or elastic network. Our result follows 
from a Landau mean-field theory for a system of confined ions in an external solution of mixed-valence 
counterions, where the ratio of squared monovalent to divalent ion concentration provides a temperature-like 
variable for the phase transition. Our analysis was motivated by long-studied volume phase transitions via ion 
exchange in ionic gels, but our findings agree with existing theory for volume-temperature phase transitions 
in charged hard-sphere models and other systems by Fisher and Levin, and McGahay and Tomozawa. Our 
mean- field model predicts a continuous line of gas-liquid- type critical points connecting a purely monovalent, 
divalent-sensitive critical point at one extreme with a divalent, monovalent-sensitive critical point at the 
other; an alternative representation of the Landau functional handles this second limit. It follows that critical 
sensitivity to ion valence is tunable to any desired valence ratio. The critical or discontinuous dependent 
variable can be the confinement volume; alternatively the internal electrical potential may be more convenient 
in applications. Our simplified conditions for ionic phase transitions to occur, together with our relatively 
simple theory to describe them, may facilitate exploration of tunable critical sensitivity in areas such as ion 
detection technology, biological switches and osmotic control. 



I. INTRODUCTION 

Phase transitions in charged networks, as manifested 
by dramatic and reversible swelling of a polymer gel, have 
been studied for over thirty years.^"^ In such gels expan- 
sion is driven by the osmotic pressure of mobile coun- 
terions that are free to exchange wi th a n external ion 
population in a Donnan equilibrium.!^^ Osmotic pres- 
sure varies under exchange of divalent for monovalent 
ions, so that that such systems can be critically sensi- 
tive to — among other things — the relative external con- 
centrations of multivalent ions. Accordingly, such phase 
transitions have long been discussed with an eye to ap- 
plications such as switches, artificial muscled and metal- 
ion detection,!^ or as candidate mechanisms for essential 
biophysical processes.!-^' As in familiar liquid- vapor tran- 
sitions, some element of self-attraction is always required 
to drive the transition; in gels, effective self-attractions 
are known to arise from relatively complicated network 
effects as theoretically described by Flory.^ In fact, such 
network behavior is by itself r ich e nough to yield critical 
behavior with no ions present.'^ESl 

Seeking a simplified framework for critical ionic sensi- 
tivity such as that observed in charged networks, we have 
recast the theory of ionic phase transitions into a solvable 
mean-field formulation, where we included in the theory 
a nonideality of Debye-Hiickel type, which is correct in 
the dilute limit for all ionic solutions. We found that 
in principle a self-attracting network is not necessary for 
a discontinuous phase transition in the presence of the 
power-law nonideality, which itself acts as an effective 
self-attraction mediating an ionically driven phase tran- 
sition of the gas-liquid type. Our model requires us to 
include a self-repulsive term to avoid runaway collapse 
due to the non-ideal term. In biological applications, in- 
tracellular charged proteins may play the role of confined 



charges, and steric exclusion among them would natu- 
rally provide the self-repulsive term. In faciei' charges on 
intracellular proteins in cells are present in densities such 
that Debye-Hiickel criticality could play a role in osmotic 
control within organisms. 

Within our mean-field theory, we find that by changing 
the magnitude of our self-repulsive term, we can move 
the system through a line of critical points. In this way, 
criticality is in principle tunable to occur at any desired 
value of external divalent ion fraction. 

Our finding of Debye-Hiickel criticality was made in 
the context of ion exchange in aqueous systems j'i^l but 
it is in accord with established work in the context of 
volume-temperature transitions in ch arged hard-sphere 
fluids initiated by Fisher and Levin,'i^IIll in semicon- 
ductor electron-hole fluids, glasses, and molten salts 
by McGahay and Tomozawa^ and in neutral polyam- 
pholytes (overall-neutral charged polymers) by Barbosa 
and Levin.li^ 

Both our simplified formulation of ionic transitions and 
our prediction of criticality in aqueous systems much sim- 
pler than gel networks may illuminate mechanisms of 
tunable critical sensitivity to ion valence or concentra- 
tion that could underlie, for example, biophysical cel- 
lular functions such as homeostasis.lJ^ The theory may 
also lend itself to engineering applications involving ion 
detection. For these purposes the internal electric poten- 
tial which we calculate, might be a more convenient 
dependent parameter than volume in applications to ion 
detection or in biophysical roles for critical ionic sensi- 
tivity. 
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FIG. 1. Schematic mechanism for a phase transition under 
changing relative concentrations of mobile monovalent and 
divalent counterions; confined charges are chosen negative in 
the diagram. Divalent ions can exchange for monovalents, re- 
ducing the internal osmotic pressure while preserving approx- 
imate neutrality. An effective self-attraction of the confined- 
charge system can lead to a volume phase transition, with 
external divalent fraction serving as an effective temperature 
parameter. 



II. DONIMAN EQUILIBRIUM WITH NONIDEALITY 



Consider a population of A'o ions each of charge go (of 
either sign, in our treatment) that may be bound to a me- 
chanical structure, such as a polymer network, or other- 
wise confined within a volume permeable to coun terion s. 
These are the conditions for Donnan equilibriurrP^E^ in 
which osmotic pressure of excess counterions goes hand- 
in-hand with an internal voltage $, relative to outside, 
with the same sign as qo- Even for $ 7^ 0, neutral- 
ity holds to a good approximation whenever the voltage 
drop $ occurs only at the boundary, or more generally if 
A^o ^ C^/qn, where C is the system capacitance. 

Osmotic pressure is exerted via the electrical poten- 
tial drop at the gel boundary. External concentrations 
of counterions may be much more dilute than the same 
species within the confinement volume. Nonetheless, 
small changes in the relative concentrations of external 
ions of differing valence can lead to large changes in in- 
ternal counterion valence ratios and thus in the internal 
osmotic pressure, as sketched in Fig. 1. When effective 
self-attraction is present in the confined-ion system, to be 
supplied in our case by the non-ideality of the internal 
solution, phase transitions in volume and a critical point 
analogous to that in a gas-liquid system can occur. 

The immobile ions, confined within a variable volume 
V, have concentration c = Nq/V = 1/v. Monovalent 
and divalent counterions with charges — —qo and 
qt = — 2go are introduced at external concentrations a 
and 6, and have within V the concentrations a' = Na/V 
and b' = N^/V. Their free diffusion in and out is con- 
trolled by the the electrical potential $. We introduce 
nonideality of mobile ions by way of a Debye-Hiickel in- 



teraction or correlation free energj^H 



3/2 



(1) 



where qa = ZaC. The essential features of F^h are that it 
is negative and contains mobile-ion concentrations raised 
to the 3/2 power. The constant Vb is proportional to the 
cube of Bjerrum length. In standard aqueous conditions. 
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To obtain closed-form solutions, we omit confined ions 
from the sum, and replace the factor (a'-t-46')^/^ by {a' + 
26')^/^, which will equal c^^^ once neutrality is imposed. 
Thus we employ the Debye-Hiickel-like interaction term 

Fint = -v'/'k^TVia' + 2b'f^\ (3) 

With a' = Na/V and h' — N^/V and taking derivatives. 
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where activity coefhcients 7a a nd for the monovalent 
and divalent ions are defined bjl™ 



^la = Ma + ksT\n{-faa). 



(7) 



Proportionality of In 7 to the square root of ionic con- 
centration is the hallmark of Debye-Hiic kel be havior, a 
theory appropriate to low ionic strength.l^^ISS! The free 
energy, with $ externally controlled, is 



F{Na,Ni,,V,^)=Fo{V) + 

F,^,{Na,Nb,V) + (iVo-TVa 
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Here ^o(^) describes mechanical constraints on the im- 
mobile ions, such as a polymer network carrying the fixed 
ions or a membrane containing them. The logarithmic 
terms are the ideal free energy of the mobile ions, and 
Co is a concentration scale which will cancel out. The 
mechanical contribution to pressure is —dF^/dV = Pg 
(equivalently a function of v or c). The overall pres- 
sure and the chemical potentials /Iq = dF/dNa and 
fib — dF/ dNb are 



P = Po{v) + {a' + b')k^T + Pint 

lad' 
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If we set fia ~ k^Tln^a/co) and /it, = k^T \n{b / cq) , a and 
b become effective concentrations, referred to an ideal 
external solution. To fix $, wc impose neutrality, c = 
a' + 2b', and introduce a dimensionless potential 



(12) 



which will always be positive, since $ will always be of 
the same sign as qq. We have altogether 

P = Pq{NqIc) + \{a' + c)k^T - \vl^'' c^'^k^T (13) 
= ln(a7a) - fi;^/^ c^'^ (14) 
2(1) = \n{b'/b) - ivl" c^/^ (15) 

which, together with c = a' + 2&' = Nq/v, we will solve to 
find P{v). Combining the equations involving 0, putting 
2b' + a' — c = for neutrality, and solving yields 



a = 



(1 + {8b/a^)c) 



1/2 



We define the dimensionless divalent parameter 
P = 
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To motivate the ratio ^ ~ by a chemical analogy, 

imagine divalent counterions B^+, complexcd with 

an A'o-valent entity 0^"+, cooperatively exchanging with 
Nq monovalent counterions A+ according to 

CA;Vo + (iVo/2)B2+ ^ CB^„/2 + iVoA+. (18) 

Equilibrium, in the sense of mass action, then gives 



Na/2 



■ (19) 



As A^o ^ oo, the monovalent-divalent exchange CB„ -H- 
CA2„ becomes discrete at a particular value of b/a^. 



III. PHASE BOUNDARY 

The dimensionless pressure p{v) = Pvs/k^T is 



P(v) = po{v) 
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where po = PoVu/k^T. To make the system stable 
against collapse to t; = (due to the interaction term) 
we require that po(^^) include a repulsion diverging faster 
than l/ii'^/^. We will use the minimal choice 



Po{v) = a (^) 
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with a dimensionless. The dimensionless electrical po- 
tential is 



(^ = ln 
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FIG. 2. Coexistence diagram in dimensionless volume and 

pressure for a system tuned by repulsion a to the monovalent 
(A = 1) critical point. The divalent parameter /3 rises from 
zero upon addition of divalent ions, so that P < is not 
possible for A = 1. In similar coexistence diagrams for < 
A < 1, the upper region will be accessible. 



Prom the Gibbs-Duhem relation (at constant T, or here 
/?) we obtain the chemical potential, linking coexisting v 
and v' at specified /?, as 



fj,{v) = pv — J p{v)dv 
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Self- intersections of the curve {p{v) , iJ,{v)) yield phase 
boundaries as in Figs. 2 and 3. 



IV. CRITICAL LINE 



We solve our model in terms of the parameter 



(24) 



At a critical point the conditions p'{v) = p"{v) = and 
p'(.t) = p" [x) = are equivalent. Regarding the critical 
value of p as a func tion of x and /?, and introducing the 
parameter A = 1/ y/l+PcxI, we have the a;-derivatives 



P:, = Xc[4.axl- Ixc + l +A] = (25) 
= 12Q!a;^-3a;c + l + A^ = (26) 
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FIG. 3. Critical line (dotted) connecting the monovalent 
(A = 1) and divalent (A = 0) critical points. The divalent 
parameter is /3 and p is dimensionless pressure. Selected crit- 
ical points and their phase boundaries (solid) are shown. The 
phase boundary becomes inaccessible (moving to j3 > oo) for 
the critical point at A = 0. 



Eliminating a in favor of A gives the line of critical points 
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The interval < A < 1 corresponds to oo > /3 > 0. Al- 
though we view a(A) as the mechanism by which ionic 
critical points can be tuned, we can also view Uc = Q!(A) 
as a critical value of the repulsion strength. With ionic 
conditions fixed, wc could drive transitions by modulat- 
ing a, analogous to transitions in ionic gels induced by 
variation of solvent composition.^ 

About each critical point, with 5v — v — Vc and (5/3 = 
l3 — jSc, we construct a Landau expansion for the pressure. 



p(/3,e) = Pc-A6p + BSPe-C( 



(31) 



where e = Sv + k5/3 is an order parameter linear in Sv 
and 6(3. In a potential application, we imagine tuning a 
system to a critical point, [3 = /3c and p = Pc- Within 
mean-field theory, when the divalent ratio is changed by 
6/3 ~ (3 — Pc one predicts a singular expansion or contrac- 
tion 



-{A/Cf/^ [6/3) 



Nl/3 



(32) 



where the cube root is taken with the same sign as 6/3. As 
discussed below, this singular behavior in volume might 



in practice be better monitored via the electric potential 
than the volume change 6v or e. 

To evaluate the Landau coeffients, we expand the pres- 
sure around the critical point, to third order in 6v = v—Vc 
and to first order in 6/3, giving 



piv,/3) w pc+p/36/3 +pyp6/3e- 



(33) 



where /3 and v subscripts denote partial derivatives and p, 
Pu, etc. are critical values. We have chosen k — pyy/3 /pwv 
to eliminate an term. Evaluating the various deriva- 
tives, the Landau parameters along the critical line are 
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where Xc and ac are known functions of A from Eqs. (28) 
C has no zeros within 



< A < 1, while K 
is zero only at A = 0. At the monovalent critical point 
iP = 0, i.e. A = 1) Eqs. ([27| through ^ yield pc = 
32/27 « 1.185, Xc = 8/3 « 2.67, ac = 9/128 « 0.0703, 
the Landau coefficients are 
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and the order parameter is e — 6v + {1/2)6/3. Here 
the system enters coexistence for any value of > 0, 
with e w — (0.158) Z?!/'^. A system tuned to the monova- 
lent critical point moves with infinite response towards 
a smaller volume with the introduction of any divalents. 
In the liquid-gas analogy, this path follows the critical 
isobar. 

At the purely divalent critical point (/3 = 00, A = 0) 
both A and B vanish, but we can here characterize the 
line of critical points by the alternative parameterization 



Pi/3,e) = pc + A'6ip-'n'B'6i/3-'n 
A' = Xcil~Xf/^/{l + X)^ 
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The parameter /3 is proportional to the monovalent 
concentration at fixed divalent concentration. The vol- 



ume singularity ( 32 ) can be rewritten as 
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where again the cube root has the same sign as its ar- 
gument. At the divalent critical point (/3 = 00, A = 0) 
Eqs. ([27| through (|30| yield Pc = 4/27 sa 0.1481, Xc = 
4/3, ac = 9/64 w 0.1406, 



A' = 4/3, B' ^ 32/27 
C = 2^^/t w 0.4162 
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For this divalent critical point at (3~'^^^ = 0, the order pa- 
rameter is simply the volume, e = Sv. A purely divalent 
system that is tuned to be critical moves with infinite 
response towards a larger volume with the introduction 
of any monovalents, as Sv « (1.214) /3^^/^. 

Returning to the dimensionless potential (j> = qo^/k^T 
in Eq. ( 22 ) , along the critical line with a constant we find 
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within a dilute screening limit, even within the confined- 
ion region. 

With regard to biological systems, we are applying the 
results here to biological mechanism^^ such as homeosta- 
sis (control of ion concentration) and cell volume control, 
and in fact find that ionic levels in cells and charges on 
intracellular proteins do fall within ranges appropriate 
for Debye-Hiickel criticality to play a role. 
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Eq. (46 ) is nonzero along the entire line of critical points. 
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V. CONCLUSIONS 

Our simplified conditions for criticality to occur in 
aqueous solutions, and our mean-field framework for an- 
alyzing them, are intended to motivate investigation of 
ionic phase transitions in ever simpler systems. As a 
fascinating analogue of liquid-gas transitions, and in the 
context of ionic gels, such transitions have been a topic 
of discussion for many years. Potential applic ation s such 
as metal- ion detection and electrical sensitivitjEEl should 
remain applicable to phase transitions as discussed here. 
Of particular interest is our finding of a line of critical 
points, such that ion-sensitive criticality is in principle 
tunable to occur at any value of the divalent fraction. 

We arrived at our finding of Debye-Hiickel criticalitjfiSI 
in the context of swelling transitions in ionic gels, a 
line of experimental and theoretical work initiated by 
T. Tanaka, but subsequently found this conclusion to 
be in accord with a distinct line of theory originating 
with Fisher and Levirfl^ for temperature-density phase 
transitions in charged hard-sphere fluids, and of McGa- 
hay and TomozawcP^ for transitions in electron-hole flu- 
ids in semiconductors, glasses, and molten salts. Sub- 
sequent theory by Barbosa and Levin^^ focused on neu- 
tral polyampholytes. Our treatment of nonideal Donnan 
equilibrium remains somewhat more phenomenological 
than that of these workers, who employed ion-pairing 
refinements of Bjerrum^"* and Fuoss^^ to obtain quan- 
titatively accurate pha se boundaries within their explic- 
itly defined models.'^ For simplicity, we have remained 
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